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ABSTRACT 

We present host-galaxy velocity dispersions of 12 local (mainly Palomar-Green) quasi-stellar 
objects (QSOs) measured directly from the stellar CO absorption features in the H band. The 
mean bulge velocity dispersion of the QSOs in our sample is 186 kms -1 with a standard deviation 
of 24 kms -1 . The measurement of the stellar velocity dispersion in QSOs enables us to place 
them on observational diagrams such as the local black-hole mass to bulge-velocity-dispersion 
relation and the fundamental plane of early-type galaxies. Concerning the former relation, these 
QSOs have higher black hole masses than most Seyfert 1 AGNs with similar velocity dispersions. 
On the fundamental plane, PG QSOs are located between the regions occupied by moderate-mass 
and giant ellipticals. The QSO bulge and black hole masses, computed from the stellar velocity 
dispersions, are of order 10 11 M and 10 8 M respectively. The Eddington efficiency of their 
black holes is on average 0.25, assuming that all of the bolometric luminosity originates from the 
active nucleus. Our data are consistent with other lines of evidence that Palomar-Green QSOs are 
related to galaxy mergers with gas-rich components and that they are formed in a manner similar 
to the most massive Ultraluminous Infrared Galaxies, regardless of their far-infrared emission. 
However, PG QSOs seem to have smaller host dispersions and different formation mechanisms 
than QSOs with supermassive black holes of 5 x 10 s — 10 9 Mq that accrete at low rates and reside 
in massive spheroids. 



Subject headings: infrared: galaxies — galaxies: active — galaxies: evolution 
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galaxies: formation 



1. Introduction 

In the current picture of galaxy formation and 
evolution, starburst and active-galactic-nucleus 
(AGN) activity are believed to be closely linked to 
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each other and to merger events (e.g., Norman & 
Scoville 1988; Haehnelt & Rees 1993; Kauffmann 
et al. 2003; Di Matteo et al. 2005; Springel et al. 
2005b; Hopkins et al. 2005; Li'pari & Terlevich 
2006). The common feature is the presence of gas, 
which is indispensible for initiating starbursts and 
feeding the AGN and which has been observed in 
local QSOs by Scoville et al. (2003) and Evans et 
al. (2001). However, the details of how mergers 
of gas-rich galaxies trigger either type of activity 
are not well understood, with the uncertainties 
mainly originating from our insufficient knowl- 
edge and treatment of the interstellar medium 
physics. Of such a nature is the debated ques- 
tion (e.g., Joseph 1999; Sanders 1999) whether 
or not quasi-stellar-object (QSO) phases can be 
associated with Ultraluminous Infrared Galaxies 
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(ULIRGs 2 ; Sanders & Mirabel 1996). 

Various scenarios have been proposed that re- 
late QSOs and ULIRGs. Sanders et al. (1998a,b) 
suggested that after the nuclei of two merging 
galaxies coalesce, the IR emission that arises from 
dust enshrouding circumnuclear starbursts and 
AGN is strong enough for the system to reach a 
QSO-like luminosity. Later in time, AGN winds 
and supernova feedback clear out the dust and gas 
from the nuclear region. The system goes through 
an optically bright phase before further accretion 
and star formation arc finally terminated. Accord- 
ing to these authors, every ULIRG should eventu- 
ally go through a QSO phase once its nucleus is 
revealed. It is, of course, possible that by the time 
the nuclear region is optically less obscured, the 
gas that remains does not fuel the AGN at a rate 
sufficient to make the latter shine as bright as a 
QSO. Therefore, a more plausible scenario is that 
some ULIRGs may evolve into QSOs, depending 
on the amount of gas consumed during their early- 
merger stages. Recent models predict a short (up 
to 10 8 yrs) QSO phase after the nuclear coales- 
cence of gas-rich galaxies (Di Matteo et al. 2005; 
Springel et al. 2005b; Hopkins et al. 2005; Cat- 
taneo et al. 2005). The outcome of these simu- 
lations depends strongly on the treatment of the 
interstellar medium. 

Another scenario is that QSOs and ULIRGs 
do not necessarily follow each other on some evo- 
lutionary sequence, but they are both triggered 
by similar conditions (e.g., Canalizo & Stockton 
2001). This assumption could hold, for example, if 
ULIRGs arc triggered by major mergers and QSOs 
by minor mergers of gas-rich galaxies (e.g., Canal- 
izo & Stockton 2001, Veilleux et al. 2006). 

All these scenarios were initially motivated 
by the fact that ULIRGs have IR luminosities 
> 10 12 Lq and number densities comparable to 
those of local QSOs (Sanders et al. 1988a) in the 
Bright Quasar Survey (BQS; Schmidt & Green 
1983) within the Palomar-Grcen (PG) catalog. 
Furthermore, PG QSOs have on average "warmer" 
IR spectral energy distributions (SEDs) than 
ULIRGs, which agrees with potential evolution- 
ary schemes (Sanders et al. 1989). Some of them 
have host galaxies with signs of recent interaction 

2 Gas-rich mergers with infrared, 8-1000 fim, luminosity 
greater than 10 Lq. 



(e.g. Surace et al. 2001). 

Canalizo & Stockton (2001) have pointed out 
that the evolutionary link between QSOs and 
ULIRGs is best elucidated by detecting starbursts 
in QSO host galaxies with the aid of spectroscopy. 
In the optical regime, spectra of the QSO host 
galaxy (or of star-forming regions in it) have been 
successfully modelled by Canalizo & Stockton 
(2001) for the derivation of the stellar ages and 
populations. Spectroscopy also provides crucial 
information about the stellar kinematics of the 
QSO host galaxy. However, this information is 
hard to extract from optical spectra, where the 
dilution of the host light by the AGN continuum 
is high. The best-suited wavelength regime for 
the extraction of the galaxy dynamics is the near- 
infrared (NIR) H band. The ratio of the host 
galaxy to the QSO photon flux is at a maximum 
there, because the SED of many stellar popula- 
tions has a maximum at ~ 1.6 ^m and the AGN 
flux has a minimum at ~ 1.2 /im (Elvis et al. 
1994). The AGN flux increases at shorter wave- 
lengths due to accretion-disk power-law emission 
and at longer wavelengths due to thermal emission 
of hot dust. Host-galaxy dynamic studies were 
performed both in the H and in the K band by 
Oliva et al. (1995) and Oliva et al. (1999). How- 
ever, these studies mainly presented observations 
of Seyfert 1 AGNs with bolometric luminosities 
lower than those of QSOs. 

We aimed to investigate the possibility that 
QSOs and ULIRGs are related by comparing the 
dynamics of the stellar populations in local PG 
QSOs to those of ULIRGs. For this purpose, 
we have carried out a European Southern Ob- 
servatory (ESO) large program 3 to acquire spec- 
troscopic data for 54 ULIRGs (including sources 
presented by Genzel et al. 2001 and Tacconi 
et al. 2002) and 12 QSOs. We have presented 
the ULIRG stellar kinematics in Dasyra et al. 
(2006a, b), hereafter Papers I and II. In this pa- 
per, we present the host dynamics of (mainly) PG 
QSOs. These Very Large Telescope (VLT) spec- 
troscopic data constitute part of a larger project 
that aims to determine whether local ULIRGs and 
QSOs can be related, through near- and mid- in- 
frared spectroscopy and imaging of their hosts; the 
program is called QUEST {Quasar and ULIRG 

3 PI Tacconi 
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Evolutionary STudy) and it is described in detail 
by Veilleux et al. (2006). 

This paper is arranged as follows. We present 
our sample and briefly summarize the observations 
and data reduction techniques in § 2. The host 
and black hole (BH) properties of the QSOs in 
our sample, as inferred from the stellar dynamics, 
are given in § 3; the dynamical results for the PG 
QSOs are then compared to those of ULIRGs in 
§ 4. We discuss the plausibility of a scenario that 
relates PG QSOs to gas-rich mergers and investi- 
gate the origin of local QSO populations in § 5. A 
summary is presented in § 6. 

2. Observations and data reduction 
2.1. Sample selection 

To investigate the evolutionary scenarios as- 
sociating ULIRGs with QSOs, we need to select 
amongst the local QSOs whose AGN properties 
most resemble those of ULIRGs. Since the AGNs 
in ULIRGs accrete at high rates (e.g., Paper II), 
the local QSOs with the most active AGNs are 
ideal candidates. The Bright QSO sample of the 
Palomar-Grecn catalog (Schmidt & Green 1983) 
is well suited to this requirement; the point- like 
appearance, the U-B cutoff selection, and, mainly, 
the B-band magnitude threshold (<16 mag) of 
these sources favors the selection of the most ac- 
tive local AGNs (Jester et al. 2005). 

The plethora of data available in the literature 
for PG QSOs allows for an optimal choice of a 
suitable sample. For example, measurements of 
the black- hole mass, Mbh, have been performed 
for 19 PG QSOs by Kaspi et al. (2000) and Pe- 
terson et al. (2004). Approximately half of the 
sources in our sample were selected to have such a 
A/bh measurement. The remaining sources were 
chosen from the Surace et al. (2001) sample of 
so-called IR-excess QSOs. These are PG QSOs 
with strong mid- and far- infrared (MIR and FIR) 
emission compared to their optical luminosity; of 
the PG sample, these are believed to most closely 
resemble ULIRGs. One more source, LBQS 0307- 
0101 (Chaffee et al. 1991), was selected from the 
Large Bright QSO Survey, since it fits well within 
the observational constraints: its £?-band magni- 
tude is similar to that of our entire sample. All 
the sources that we observed were at redshift z < 
0.1. This upper redshift cutoff was selected to 



be smaller than that of the Spitzer subsample of 
QUEST, ^0.3, to maximize the host signal and 
facilitate the extraction of the stellar kinemat- 
ics. The identifiers, coordinatess, redshifts, and 
B-band magnitudes of the 12 sources selected and 
observed are given in Table 1. 

To investigate how well our PG subsample rep- 
resents the entire BQS PG population, we ex- 
amined the optical and MIR/FIR luminosities of 
our sources as tabulated in Table 2. The opti- 
cal luminosities XL\ (5100^1) are from Peterson et 
al. (2004) and Vestergaard & Peterson (2006). 
We computed the infrared luminosities Ljr from 
the 12, 25, 60, and 100 /zm fluxes (Sanders et al. 
1989; Moshir et al. 1990; Haas et al. 2003) us- 
ing the Sanders & Mirabel (1996) prescription. 
The IR luminosity of all sources in our sample is 
greater than 10 11 L Q . This luminosity threshold 
defines Luminous Infrared Galaxies (LIRGs). In 
Table 2, we also present the ratio of the IR to op- 
tical/ultraviolet luminosity of our sources, as de- 
rived from Surace et al. (2001) and Guyon et al. 
(2006). The latter is denoted as "big blue bump" 
luminosity, Lbbb, and it is computed by integrat- 
ing the luminosity per unit wavelength from ~0.01 
to 1 /im. 

In Fig. 1, left panel, we plot the optical ver- 
sus the infrared luminosity of PG QSOs. To con- 
struct this diagram, we used PG QSOs in the en- 
tire QUEST redshift range. We also used those 
PG QSOs of the Peterson et al. (2004) and Vester- 
gaard & Peterson (2006) samples that have a flux 
measurement, instead of an upper limit, at 12, 25, 
60, and 100 /im. The sources with a mid- or far- 
infrared flux detection, which comes either from 
Infrared Astronomical Satellite (IRAS) or from In- 
frared Space Observatory (ISO) data, are ~ 50% 
of the sources with with a \L\(5100A) estimate. 
The PG QSOs of our sample are plotted as filled 
circles and all others as open circles. Both the 
optical and IR luminosities of the sources in our 
sample are very close to those of the mean of the 
IR-detected PG QSO population, indicating that 
our sources are representative of the most common 
PG QSOs but do not span the full range of AGN 
luminosities. In Fig. 1, right panel, we plot the 
60-25 /im versus 100-60 /im color index a diagram 
(Dc Grijp et al. 1987) of local PG QSOs: the sym- 
bols used are identical to those in the left panel of 
Fig. 1. The MIR color indices of our sources are 
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typical of all QSOs used in this diagram, again 
indicating that the sample we compiled is repre- 
sentative of the local IR-detected PG QSOs. 

2.2. Data acquisition and analysis 

Our long-slit spectra were obtained with the 
ISAAC spectrometer (Moorwood ct al. 1998) 
mounted on the Antu telescope unit of the ESO 
Very Large Telescope on Cerro Paranal, Chile. 
The QSO observations were performed in both 
service and visitor mode typically under excellent 
seeing conditions (optical seeing of ^0'.'5-0'.'6) to 
minimize the effects of AGN light diluting the host 
signal. The integration time for each exposure was 
300 s to avoid saturation of the detector, and the 
total on-source integration time varied from 160 to 
320 mins. The slit width was Of! 6 and the detector 
scale is 0'.'146 per pixel (see Paper I). The central 
wavelength was in the H— band and varied from 
1.70 to 1.77 /im, depending on the redshift of each 
source and the wavebands of high atmospheric 
transmission. The spectral resolution R = A/AA 
was 5100. The slit position angles (PAs) and the 
respective integration time per PA are tabulated 
in Table 1 for the 12 sources observed. The red- 
shift range of these QSOs (0.050 < z < 0.112) is 
on average lower than that of the ULIRGs in our 
sample (0.018 < z < 0.268), and the integration 
times ~2 times longer (on average), to ensure a 
high signal-to-noise (S/N) ratio of the host-galaxy 
signal. 

For most of the QSOs in this study, the 
H— band effective radius i? e ff and mean surface 
brightness n c g within it are measured from Hub- 
ble Space Telescope (HST) NICMOS imaging 
data (McLeod & Rieke 1994a; Veilleux et al. 
2006). For some sources, NIR imaging data (as- 
sisted by adaptive optics) are also available from 
Guyon et al. (2006) 4 . The quantities i? e ff and 
/U e fi are given in Table 1. The R c ff values are 
converted to the cosmology used in this paper 
(H a =70 km s- 1 Mpc" 1 , O ro =0.3, fi t otai=l). 

The data reduction tasks and procedures used 
to extract the host-galaxy spectra from the ac- 
quired ISAAC data are described in detail in Pa- 

4 It is not possible to use our acquisition images to derive 
R e ff> as we did in Papers I and II for the ULIRGs, since 
the photon counts in our short-exposure-time acquisition 
images are strongly dominated by the point- like AGN, lead- 
ing to i? c ff values biased towards small radial extents. 



pers I and II. In the rest of this subsection, we 
mainly list the differences in the reduction method 
and the derived results between the ULIRG and 
the QSO data. 

The main difference originates from the fact 
that photon counts associated with the AGN con- 
tinuum are roughly one order of magnitude greater 
than those of the host galaxy, leading to a signifi- 
cant suppression of the stellar absorption features. 
For this reason, the highly nucleated AGN emis- 
sion needs to be avoided for the accurate extrac- 
tion of the stellar kinematics. To optimize the host 
S/N ratio per slit, we extracted the spectra from 
the widest possible aperture (typically ±l'.'0-l'/4 
from the center), excluding the very central re- 
gion (±0'.'3-0'/4). We then combined the results of 
the two slits into a single spectrum per object. 

The method of determining the stellar velocity 
dispersion a from the reduced spectra was based 
on a Fourier correlation quotient technique (Ben- 
der 1990). This method uses a template spectrum 
to provide the Dopplcr-broadened profile of the 
stellar kinematics along the line-of sight (LOS). 
We fitted a combination of a Gaussian and a low- 
order polynomial to the broadened profile to deter- 
mine the projected stellar velocity dispersion. The 
fit was performed to each bandhcad individually 
and the quoted errors are the standard deviation 
of all measurements. In contrast to Papers I and II 
where a was measured from the central aperture, 
the stellar dispersion derived for the QSO hosts 
was typically measured at a radius of ^i? c ff/4. In 
simulations of gas-rich mergers, the velocity dis- 
persion at the center and at the effective radius 
of the remnants tend to differ by ~10% (Bendo & 
Barnes 2000). We ignored aperture effects since 
the expected differences are typically within our 
error bars. We could not extract information on 
the rotation of the host galaxy for the QSOs, given 
that it was necessary to perform pixel rebinning to 
achieve the maximum possible S/N ratio. 

We observed four template stars for the extrac- 
tion of the host velocity dispersion: an M0III giant 
(HD 25472), an MHab supergiant (HD 99817), a 
K5Ib supergiant (HD 200576), and a KOIab super- 
giant (HD 179323). The spectra of the first three 
stars were presented in Genzel et al. (2001). The 
KOI supergiant was observed with the integral- 
field-unit SINFONI on the VLT, which has a 
spectral resolution i?=3000 in the iJ-band. The 
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2-d images were flat-ficldcd and converted into 
a wavelength-calibrated 3-d data cube using the 
SINFONI data cube reduction package "SPRED" 
(Schreiber et al. 2004). After the extraction of 
the 1-d spectrum from the 3-d data cube, the stel- 
lar spectrum was sky-background subtracted and 
telluric features were accounted for by compari- 
son with a stellar spectrum of known type. The 
methodology followed to apply these corrections to 
the 1-d spectrum was identical to that applied for 
the ULIRGs in our program. The ii-band spec- 
trum of HD 179323 is presented in Fig. 2. Un- 
like most ULIRGs, which are best described by 
the MOIII giant (see Papers I and II), some PG 
QSOs are best fit by the KOI supergiant. In some 
cases, a linear combination of the two provides the 
best description of the QSO-host spectrum. The 
template used for each source and the measured 
velocity dispersion are given in Table 3. The use 
of different stellar templates leads to a measure- 
ments that differ typically by 5 to 20%. 

We plot the restframe spectrum of each QSO 
host galaxy in Fig. 3. For each source, we over- 
plot the selected stellar template after convolving 
it with a Gaussian of dispersion equal to that mea- 
sured. The spectra shown in Fig. 3 are convolved 
to the resolution of the SINFONI observations for 
the sources that are best described by the KOI su- 
pergiant or the combination of the KOI and the 
MOIII templates. 

The velocity dispersion measurements can be 
affected by the strong dilution of the host-galaxy 
spectrum by the AGN continuum, which alters the 
shape of the wings of the stellar absorption fea- 
tures. To quantify this effect, we have constructed 
artificial QSO ISAAC data from which we mea- 
sured the bulge velocity dispersion as we did for 
the real data. We convolved the spectrum of our 
MOIII stellar template with a Gaussian of a =200 
kms -1 . We scaled its photon counts to those of 
a de Vaucouleurs bulge with the mean iJ-band 
magnitude and the mean redshift of the sources 
in our sample, to#=13.5 mag (McLeod & Rieke 
1994a; Veilleux et al. 2006) and z = 0.076 respec- 
tively, that falls in a 0'.'6 slit under seeing condi- 
tions of 0'.'6. We diluted this host-galaxy signal 
with a continuum whose strength was determined 
by counting the photons that originated from a 
point source of (mean) m# = 13.4 mag and that fell 
in the slit under the assumed seeing conditions. 



After adding Poissonian noise, we extracted the 
bulge dispersion in the same way we did for the 
real data. We repeated this procedure for 10 iter- 
ations and measured a mean a of 179 kms -1 with 
a standard deviation of 21 kms -1 . Therefore, ve- 
locity dispersion underestimates of order 10% due 
to the strong dilution of the host absorption fea- 
tures would not be unexpected. In addition to 
this, any differences between the stellar spectrum 
used to correct for the telluric features and the 
actual atmospheric absorption lines could also af- 
fect our spectra (see Paper I). This could happen 
if, for example, the residuals from the removal of 
the stellar features are at ~0.1% continuum level, 
since the QSO-host-galaxy emission is only ~ 10% 
of the AGN continuum. Both these effects would 
lead to an artificial mismatch between the tem- 
plate and the galaxy spectrum, which is effectively 
included in the uncertainties of the measured ve- 
locity dispersion. In some cases, they lead to a 
error bars that are > 50 kms -1 . 

The S/N ratio was insufficient for the detec- 
tion of stellar absorption from the host galaxy of 
PG 1211+143, probably due to a very strong di- 
lution of the host-galaxy light by the AGN con- 
tinuum. PG 1426+015 is an interacting galaxy 
with two components. The data for both the 
bright, north-east (NE) nucleus and the secondary, 
south-west (SW) nucleus are presented here (see 
Figs. 3, 4). The H- band spectrum of the NE 
nucleus is mainly dominated by the AGN contin- 
uum and it is ~30 times brighter than that of the 
fainter nucleus. 

3. Results 

3.1. Host-galaxy dynamical properties 

The stellar velocity dispersions of all QSO host 
galaxies in our sample are presented in Table 3. 
The dispersion of the SW nucleus of PG 1426+015 
(which is not a QSO) is excluded from all statis- 
tical analyses of this paper. The mean value of 
a is 186 kms -1 , with a standard deviation of 24 
kms -1 and a standard error (uncertainty of the 
mean) of 7 kms -1 . Most of the sources of Ta- 
ble 3 with a < 200 kms -1 show spiral structure or 
tidal tails (PG 0050+124, PG 1119+120, PG 1126- 
041, PG 1229+204, PG1426+015, PG 2130+099) 
in their NIR images (Surace et al. 2001; Veilleux 
et al. 2006). 
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For a stellar system characterized by a disper- 
sion velocity a, the bulge mass is computed from 

to = 1.40 x 10 6 cr 2 i? c ff , (1) 

where a is in units of kms -1 , R e g is in kpc and to 
is in Mq (see Tacconi et al. 2002; Paper II). Us- 
ing this formula we find that the mean bulge mass 
of all QSOs in our sample is 2.09 x 10 11 M©. The 
bulge mass is somewhat lower than the dynami- 
cal mass, which also takes into account the stellar 
rotation. However, the fact that the mass scales 
with 3 x (T 2 and only V TOt 2 (Tacconi et al. 2002) 
and that any disk structure observed in the QSO 
hosts is not as prominent as the bulge (Veilleux et 
al. 2006) is a good indication that the dynami- 
cal mass will be similar to that of the bulge. This 
result implies that the hosts of PG QSOs are typ- 
ically 1.5to* galaxies, for to* = 1.4 x 10 11 M 
(see Genzel et al. 2001, Paper II, and references 
therein). 

In Fig. 5 we place our QSOs on the if -band fun- 
damental plane (FP) of early-type galaxies (Djor- 
govski & Davis 1987; Dressier et al. 1987) using 
the host effective radii and mean surface bright- 
nesses from Veilleux et al. (2006) and McLeod & 
Rickc (1994a) tabulated in Table 1. We converted 
H- into K- band host magnitudes using the Two 
Micron All Sky Survey (2MASS) mean H- K=Q2 
color index correction (Jarrett et al. 2003) 5 . The 
data for the early-type galaxies are taken from 
Bender et al. (1992), Faber et al. (1997), and 
Pahre (1999) and converted to our cosmology. 
Giant boxy ellipticals (squares) are located on the 
upper-right part of the fundamental plane R e g- 
a projection, while moderate-mass ellipticals (cir- 
cles) occupy the central and the lower-left parts 
(see Fig. 5, left panel). PG QSO hosts (stars) lie 
between moderate-mass and giant Es, but closer to 
the former, both in the a-R c s projection that re- 
lates the dynamical properties of the systems (up- 
per left panel) and in the 3-dimensional view of the 
plane that takes into account their photometric 
properties (middle and lower panels). On the i? e ff- 
/i ff projection (upper right panel), also known 
as the Kormendy relation (Kormendy 1985), PG 

5 We do not discriminate between elliptical or spiral hosts 
since the difference is small for the 2MASS galaxies and 
since the hosts of some of the QSOs in our sample show 
patterns of spiral structure (e.g. Surace et al. 2001; Guyon 
et al. 2006). 



QSO hosts have only a small overlap with QSOs 
that are hosted by giant Es (Dunlop et al. 2003; 
see § 5.2). 

3.2. Black hole properties 

We compute the black hole masses of the QSOs 
in our sample by using the tight correlation be- 
tween black hole mass and bulge velocity dis- 
persion, the A'/bh — cr relation (Gebhardt et al. 
2000a; Ferrarese & Merritt 2000). We use the 
Tremaine et al. (2002) formula, 

M BH = 1.35 x 10 8 [cr/200] 4 - 02 M , (2) 

as in Papers I and II. The application of the 
A^bh — cr relation yields an average black hole 
mass of 1.12 x 10 8 M© (see Table 3). 

While the Mbh — c relation was established us- 
ing mostly local quiescent galaxies, it was subse- 
quently shown that a similar relation exists in low- 
luminosity AGNs, namely type 1 Scyfcrts (Geb- 
hardt et al. 2000b; Ferrarese et al. 2001; Nelson 
et al. 2004; Onkcn et al. 2004). The AGN black 
hole masses were measured from reverberation- 
mapping experiments (Blandford & McKee 1982; 
Peterson 1993; Peterson et al. 1998; Kaspi et al. 
2000; Peterson et al. 2004; Kaspi et al. 2005; 
Vestergaard & Peterson 2006). Reverberation 
mapping measures the size of the broad-line region 
(BLR) from the light-travel time delay between 
continuum and emission-line flux variations. Un- 
der the assumption that the BLR gas is virializcd, 
the central black hole mass is given by 

M BH - /i?BL R AU 2 /G, (3) 

where -Rblr is the size of the BLR as measured 
by the time delay, AV is the emission-line width, 
G is the gravitational constant, and / is a dimen- 
sionless scaling factor that encapsulates uncertain- 
ties in the structure and inclination of the BLR. 
The mean value (/} = 5.5 (Onken et al. 2004), 
is statistically determined by assuming that the 
A/bh — f relation is the same in both active and 
quiescent galaxies. According to Vestergaard & 
Peterson (2006), use of this single value for / 
leaves residual scatter around the Mbh — a rela- 
tion that indicates that the reverberation-based 
black hole masses are accurate to a factor of ~ 3. 
However, in individual cases the actual value of / 
depends on currently unknown factors, such as the 
inclination of the BLR (Collin et al. 2006). 
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Reverberation BH masses are available in the 
literature (Peterson et al. 2004) for four of the PG 
QSOs in our sample that have a host galaxy ve- 
locity dispersion measurement (see Table 3). We 
place these sources on the AGN Mbh — o relation 
in Fig. 6, where the QSOs (indicated by stars) are 
plotted over the AGNs (circles) of Onken et al. 
(2004) and Nelson et al. (2004) 6 . To avoid uncer- 
tainties in the mean value of /, we use the virial 
product Mbh/(/) of Eq. 3. The solid and dashed 
lines correspond respectively to the Tremaine et al. 
(2002) and Ferrarese (2002) fits, scaled down by a 
factor of (/) to match the AGN datapoints (Onken 
et al. 2004). We find that one of our sources, 
PG 1229+204, falls on the AGN relation, show- 
ing good agreement between the dynamically- 
determined and the reverberation-based black hole 
mass measurements. However, three of the four 
reverberation-mapped QSOs lie above the locus of 
Seyfert 1 AGNs with similar velocity dispersions. 
Thus, the average dynamically-determined black- 
hole masses for PG QSOs tend to be smaller than 
those measured by reverberation experiments. 
This result might be attributable to incorrect as- 
sumptions, measurement errors, or simply small 
number statistics. We consider each of these in 
turn. 

By incorrect assumptions, we are referring prin- 
cipally to the underlying assumption that a single, 
statistically determined value of / is equally good 
for all AGNs. The apparent differences can also 
reflect real differences in the host and black-hole 
properties between lower-mass and higher-mass 
black holes, such as those claimed for quiescent 
galaxies by Lauer et al. (2006). Alternatively, the 
offset may be indicating differences due to inclina- 
tion effects (Collin et al. 2006) or to a currently 
misestimated calibration of the factor /. This pos- 
sibility is implied by Fig. 6: the virial products 
of the PG QSOs in our sample lead to an offset 
above the AGN Mbh — o relation that brings the 
sources on the relation for quiescent galaxies with- 
out the use of a factor / (in other words, with 
(/) ~ !)• O ur result implies that the mean factor 
/ could be smaller than 5.5. 

PG QSOs could also deviate from the AGN 
Mbh — c relation due to measurement errors ei- 



For the sources with two <r measurements we use the aver- 
age of the two. 



ther in their stellar velocity dispersions or in their 
M BH estimates. We suspect that PG 2130+099 
represents a case of the latter. The reverberation- 
based mass significantly exceeds the mass pre- 
dicted by Vestergaard & Peterson (2006) for the 
optical luminosity and H(3 flux of this source (see 
below). The object is also an outlier in the BLR 
radius — luminosity relation (Kaspi et al. 2000; 
Bentz et al. 2006); the BLR size, based on rever- 
beration measurements, appears to be much too 
large for its luminosity. An overestimated value of 
the BLR radius is plausible for this object since 
its light curve was not particularly well sampled 
(Kaspi et al. 2000). 

It is also possible that there are systematic er- 
rors in the measurement of a, e.g., from dilu- 
tion of the stellar light by the AGN continuum. 
However, such differences can typically be of or- 
der ~10% (see § 2.2). By themselves, they seem 
unable to account for the discrepancy between 
the dynamically-determined and the reverberation 
A#bh estimates, which on average is a factor of 7. 
However, other systematics, such as stellar popu- 
lation effects might also contribute to the discrep- 
ancy. Errors of this type could be introduced by 
the use of NIR velocity dispersion measurements 
for the construction of the AGN Mbh — o - rela- 
tion, which is based on optical a values. Discrep- 
ancies of order 10% have been reported in the lit- 
erature between CO and various optical dispersion 
measurements, mostly for quiescent galaxies that 
contain disky structures (e.g., Silge & Gebhardt 
2003; Oliva et al. 1999). However, the corrections 
implied from the literature are not only small, but 
also have conflicting signs. It is very likely that 
they reflect differences in the methodology, spec- 
tral resolution, and stellar templates that have 
been used by the various authors. At this point, 
the possibility that there are true systematic dif- 
ferences between a measurements in the optical 
and NIR cannot be discounted. Further data are 
required to either confirm or quantify such effects. 

A third possible explanation is that the appar- 
ent displacement of the PG QSOs from the AGN 
Mbh — & relation is simply an artifact of small- 
number statistics. As a test of this possibility, we 
use a secondary method to estimate the BH masses 
for 10 objects in our sample. Reverberation- 
mapping results show a strong correlation between 
BLR radius and the luminosity of the AGN (Kaspi 
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et al. 2000; 2005; Bentz et al. 2006); it is therefore 
possible to infer the BLR radius and the black- 
hole mass of the AGN by making use of an optical 
(or UV) flux measurement and a "single-epoch" 
spectrum that provides emission-line widths (e.g., 
Vestergaard 2002). Using such relations, Vester- 
gaard & Peterson (2006) have recently provided 
Mbh estimates for PG QSOs without direct re- 
verberation measurements. Of those, we use the 
values computed from the width of the H/3 line and 
the continuum luminosity at 5100 A. The statisti- 
cal error bars 7 that we ascribe to these black-hole 
mass estimates are based on the 1-sigma scatter 
around this luminosity-black hole relation, which 
amounts to 0.52 dex or a factor of 3.3 (Vester- 
gaard & Peterson 2006). In Fig. 7, we plot the 
10 PG QSOs (shown as stars) with single-epoch 
Mbh estimates over the Seyfert 1 AGNs of Onken 
et al. (2004) and Nelson ct al. (2004) (circles). 
We also overplot the M BH ~ 10 6 M Q AGNs of 
Barth et al. (2005) (triangles) to simultaneously 
display the behavior of the Mbh — c relation at 
the high and low-mass ends. For all sources, we 
use the Mbh/C/} estimates with (/) = 5.5, as in 
Fig. 6. By increasing the number of sources on the 
AGN Mbh — o diagram, we find a better agree- 
ment between the Seyfert 1 AGNs and the PG 
QSOs. The scatter in the relation seems to be 
larger at both the low and high mass ends than 
for the reverberation-mapped AGNs, due to the 
additional scatter that is intrinsic to the single- 
epoch scaling relations. 

On the basis of Fig. 7, it appears that the re- 
sult found for the reverberation-mapped PG QSOs 
is at least partially ascribablc to small number 
statistics, and that any correction to (/) is rather 
small. Moreover, the offsets between the single- 
epoch and the dynamically-determined values of 
Mbh are comparable both at the high and low 
mass ends of the Mbh — o relation; specifically, 
< 51og[M BH ]> is 0.15 and -0.25 for the PG QSOs 
and the low-luminosity Barth et al. (2005) sam- 
ples, respectively. That the offset is of the same 

7 Vestergaard & Peterson (2006) provided statistical uncer- 
tainties for this relation at various confidence levels for the 
sample as a whole and measurement-related error bars for 
individual sources. The statistical uncertainties are signifi- 
cantly higher than those of the flux measurements, leading 
to actual errors in the individual Mbh estimates that can 
be as high as an order of magnitude (Kelly & Bechtold 
2006). 



sign in both ends of the relation can be traceable 
to the fact that some of the QSOs in our sample 
and those of Barth et al. (2005) are narrow-line 
Seyfert 1 (NLS1) galaxies, whose black holes may 
be related in a different way to the host galaxy 
(Grupe & Mathur 2004; Barth ct al. 2005). That 
the offset between these two samples and that of 
the Onkcn ct al. (2004) AGNs is small is an in- 
dication that the relation that connects the bulge 
dispersion to the black-hole mass is global, apply- 
ing to AGNs with M B h from -10 5 to ~ 5x 10 8 
Mq. However, the measurement of a in larger 
samples of reverberation-mapped AGNs is a cru- 
cial step towards understanding where differences 
may come from and presenting a new best fit to 
the AGN data. 

To estimate at what rates our sources accrete, 
we calculate their Eddington black hole masses, 
Ms#(Edd), from their Eddington luminosities, 

-^Edd) as 

L E dd/L Q = 3.8 x lO 4 M BH (Edd)/M . (4) 

We consider that all of the QSO bolomctric lumi- 
nosity comes from the active nucleus (in contrast 
to Paper I and Paper II where we assigned half of 
the infrared luminosity of ULIRGs to ^Edd)- We 
compute the bolomctric luminosity from the opti- 
cal luminosity (see Table 1) as 

L bol = C*AL A (5100i), (5) 

where the conversion factor C is — 9 (Kaspi et 
al. 2000). The value of C can significantly dif- 
fer from one QSO to another; therefore, the accu- 
racy of its mean is rather small (within a factor 
of 3; Sanders et al. 1989; Elvis et al. 1994; Net- 
zer 2003). The estimated Eddington accretion 
efficiencies, r?Edd = MBH(Edd)/MBH, carry this 
uncertainty. Furthermore, since part of the emis- 
sion originates from the host, the values of r^Edd 
given in Table 3 should be treated as upper limits. 
The mean Eddington accretion efficiency of the 
QSOs in our sample is 0.25, in good agreement 
with McLeod et al. (1999). 

4. ULIRG vs. (IR-excess) PG QSO dy- 
namics 

To properly investigate whether QSOs and 
ULIRGs may be related, it is of particular im- 
portance to identify those QSOs with IR-excess 
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emission and compare the dynamics of the latter 
to those of ULIRGs. Several criteria have been 
used in the literature to identify IR-excess QSOs 
(see Sanders et al. 1989, Surace et al. 2001, Lipari 
1994, Canalizo & Stockton 2001, and references 
therein). Some were based on the 60 /xm flux and 
the 25/60 /im flux ratio. Sanders et al. (1989) and 
Surace et al. (2001) defined IR-excess QSOs to be 
those sources with Lm/Lbb > 0.46. Lipari (1994) 
and Canalizo & Stockton (2001), among others, 
used the diagram of MIR color indices (Fig. 1, 
right panel) to select sources in a possible tran- 
sition phase between starburst and QSO activity. 
The classification is based on the fact that AGN- 
dominated (power law) sources occupy a different 
locus than starbursts, whose characteristic inter- 
stellar dust temperature is 10 — 10 2 K (Downes 
& Solomon 1998). We opt to use the definitions 
that are based on all IRAS bands, since they are 
less susceptible to the thresholds set, which arc 
sometimes subjective and specific to the sample 
under examination. We refine the IR-excess def- 
inition to take into account both the ratio of the 
IR over the optical/ultraviolet luminosity and the 
MIR color-index diagram. As IR-excess we clas- 
sify those sources that have Lm/ibb > 0.46 and 
lie between the AGN- and starburst- dominated 
loci on the diagram of MIR color indices. 

Five sources satisfy our set of IR-excess criteria, 
namely, PG 0050+124, PG 1119+120, PG 1126- 
041, PG 1426+015, and PG2130+099. The mean 
velocity dispersion of this sample is 180 kms -1 
with a standard deviation of 13 kms -1 and error 
of 6 krns" 1 . Statistically, the difference between 
this value of a and the mean of our entire sample 
is smaller than 1 Gaussian sigma, so we find no 
particular reason to separately study the dynami- 
cal properties of IR-excess QSOs and all PG QSOs 
in our sample. The comparison that follows refers 
to our sample as a whole. 

The mean velocity dispersion of the QSO hosts 
is 1.16 times that of the ULIRG remnants in our 
sample, 161 kms _1 (Paper II). This difference in 
the dispersion can account for differences in the 
masses of ULIRGs and QSO hosts of the order 
~5x 10 Mq. However, the observed difference 
is ~ 10 11 Mq, since the mass of the ULIRG rem- 
nants is 8.91 xlO 10 Mq (Paper II) and that of the 
QSO hosts is 2.09 xlO 11 M Q . That the measured 
difference is that large can mainly be attributed 



to the effective radii of QSOs which are roughly 
twice as large as those of ULIRGs; the mean value 
of R s is 3.9 kpc and 2.2 kpc for the QSOs and 
the ULIRG remnants in our sample respectively 8 . 

We place both PG QSOs and ULIRGs on 
the fundamental plane of early-type galaxies 
(Fig. 8), to illustrate the similarities and differ- 
ences amongst their galaxy properties. To better 
visualize the region that each population occu- 
pies, we construct various 2-d and 3-d views of 
the plane, as in Fig. 5. Local ULIRGs (Gen- 
zel et al. 2001; Tacconi et al. 2002; Paper II; 
Rothberg & Joseph 2006) are plotted as trian- 
gles and PG QSOs as stars. For comparison, 
we also use local LIRGs (Shier & Fischer 1998; 
James et al. 1999; Rothberg & Joseph 2006; 
Hinz & Rieke 2006), which are plotted as circles. 
Other, morphologically-selected merger remnants 
that are not IR-bright (Rothberg & Joseph 2006) 
are plotted as boxes. The latter category can in- 
clude mergers of various types (e.g. spiral-spiral, 
elliptical-spiral, elliptical-elliptical) ; it can also in- 
clude former ULIRGs and LIRGs whose starbursts 
have now faded away. For viewing clarity we do 
not overlay any early-type galaxies on this figure. 
The position of PG QSOs on the dynamical a-R c s 
projection shows that they lie at the upper end of 
the locus of ULIRGs and LIRGs and in the middle 
of the locus of morphologically-selected mergers; 
therefore a subset of all (U)LIRGs and PG QSOs 
have identical dynamical properties. 

5. Discussion 

5.1. Relation of PG QSOs to gas-rich 
mergers 

In this section, we place our results on the host 
kinematics of PG QSOs into a wider evolutionary 
framework by combining multiwavclcngth spectro- 
scopic and imaging evidence that links PG QSOs 
to mergers or interactions of galaxies that possess 
gas. Some sources show clear signs of interaction, 

According to Tacconi et al. (2002), the measured effective 
radii of the gas-rich merger remnants may increase when 
the ultraluminous phases cease, since the strongly nucle- 
ated starburst emission will eventually fade away. Under 
this hypothesis, the difference in the stellar mass of the 
QSO hosts and the elliptical galaxies that ULIRGs will fi- 
nally form may be smaller than that measured at present 
time. 
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such as disturbed morphology (e.g. PG 1211+143; 
PG 2214+139; Surace et al. 2001; Guyon et al. 
2006) and tidal tails (e.g. PG 0007+106; Surace 
et al. 2001). Faint blobs, which could be iden- 
tified as secondary nuclei, appear in the images 
of PG 2130+099 and PG 0007+106 (Veilleux et 
al. 2006; Guyon et al. 2006. Various imaging 
analyses indicate that several sources in our sam- 
ple are unambiguously binary systems. Surace et 
al. (2001) found an elongated companion north- 
west of the main nucleus of PG 1119+120. Ac- 
cording to Veilleux et al. (2006) and Guyon et 
al. (2006), PG 1126-041 also interacts with a 
galaxy at a projected distance of 6.6 kpc. Stock- 
ton (1982) and Canalizo & Stockton (2001) 
spectroscopically confirmed that PG 0050+124 
has a companion galaxy 16.5 kpc to the west of 
the QSO. Hubble Space Telescope (HST) imag- 
ing data presented by Schade et al. (2000) indi- 
cate the presence of a second nucleus south-west 
of the bright nucleus of PG 1426+015 (also see 
Fig. 4). Our absorption- line spectroscopy of the 
host of PG 1426+015 confirms that the redshift of 
the second nucleus is identical to that of the QSO. 
The projected nuclear separation of the interact- 
ing galaxies is 4.4 kpc. PG 1426+015, PG 1126- 
041, and PG 0050+124 are amongst the objects 
in our sample for which the QSO phase has been 
reached already before the individual nuclei coa- 
lesce 9 (see also McLeod & Rieke 1994b). There- 
fore, QSO phases are not necessarily linked to the 
end of the merger or interaction process. 

Various pieces of evidence for star-formation 
exist in several of the PG QSOs in our sam- 
ple. Canalizo & Stockton (2001) have modeled 
the optical spectra of 9 ULIRGs and QSOs sus- 
pected to be in a possible ULIRG/QSO transi- 
tion phase (from their MIR color-index diagram) . 
They find that the spectra of most of their ob- 
jects are described by a combination of old and 
recently formed stellar populations. The ages de- 
rived for the starbursts were < 300 Myr. Ac- 
cording to Canalizo & Stockton (2001), one 
of the sources in our sample, PG 0050+124 (or 
I Zw 1) has ongoing star formation (Schinnerer et 
al. 1998). In the MIR, Spitzer IRS spectroscopy 
of QSOs and ULIRGs in our QUEST project (PI 

'On the other hand, PG 1119+120 does not officially fulfill 
the QSO definition, since its bolometric luminosity is < 
10 12 L Q (Surace et al. 2001). 



Veilleux) indicates the presence of starbursts in 
PG QSOs. Polycyclic Aromatic Hydrocarbon 
(PAH) molecule emission has been individually 
detected for 11 of the 26 QSOs observed; for the 
remaining sources, stacking of the spectra also 
reveals the presence of PAHs (Schweitzer et al. 
2006). Based on the strength of the PAH emission, 
Schweitzer et al. (2006) find that star formation 
in PG QSOs is much stronger than what indicated 
by emission lines at shorter wavelengths, such as 
the [Oil] line at 3727 A (Ho 2005), which can be 
affected by obscuration. 

The presence of gas, which is required to initi- 
ate star-formation and AGN-accretion events, has 
been confirmed in several PG QSOs. Molecular 
CO and HCN J=(l-0) emission lines have been de- 
tected by Evans et al. (2001; 2006), and Scoville et 
al. (2003) in a total sample of 18 PG QSOs. The 
CO luminosities were at most an order of mag- 
nitude lower than those of ULIRGs (Evans et al. 
2001; 06; Downes & Solomon 1998 and references 
therein). Converting the CO luminosity into gas 
mass implied that PG QSOs can have molecular 
gas reservoirs as massive as ^10 10 M Q . For this 
computation, a Galactic CO-to-H2 conversion fac- 
tor of 4 MQ l K^ 1 [km/ 's]~ 1 pc~ 2 was used. How- 
ever, this gas mass may be an upper limit, if the 
conversion factor for PG QSOs is closer to that 
of ULIRGs, 0.8 MeK-^km/s^pc- 2 (Downes & 
Solomon 1998), than to that of the Milky Way, 
due to a potentially compact, high-pressure gas 
distribution in the QSOs. 

In some cases, the evidence that relates lo- 
cal PG QSOs and ULIRGs is strong. For ex- 
ample, a few sources are known to fulfill both 
the ULIRG and QSO classification criteria (e.g. 
PG 0157+001; PG 1226+023). More often, the 
IR luminosity output is similar to that of LIRGs. 
This result is still consistent with an evolutionary 
scenario suggesting that one population evolves 
into another, since some LIRGs are plausibly the 
progenitors or descendants of ULIRGs. However, 
the LIRG population as a whole is believed to orig- 
inate from a wider variety of interactions or merg- 
ers, e.g., from mergers of larger progenitor mass 
ratios than those of ULIRGs (Ishida et al. 2007, 
in preparation). This is an argument against a 
scenario where all ULIRGs undergo a QSO phase 
(and vice- versa), and in favor of a scenario where 
QSO activity is triggered by progenitors or ini- 
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tial conditions that are similar but not identical 
to those of ULIRGs. Examples include some PG 
QSOs in our sample (PG 1126-041, PG 1229+204, 
and PG 2130+099) that are amongst the sources 
with the lowest optical luminosity and that have 
signatures of spiral structure (Surace et al. 2001; 
Vcillcux et al. 2006). If these sources have indeed 
gone through a recent interaction, they may have 
been involved in unequal-mass major (^3:1) or mi- 
nor (>3:1) mergers, since the latter are known to 
form remnants with significant angular momen- 
tum (e.g., Bendo & Barnes 2000). Therefore, it 
seems that both scenarios about the relation of 
ULIRGs and QSOs are partially valid. 

5.2. Relation of PG QSOs to other QSO 
populations 

It is important to point out that local gas-rich 
mergers are not the direct progenitors of all local 
QSOs, and that the PG QSO sample does not in- 
clude objects representative of all local QSO popu- 
lations (Vestergaard & Peterson 2006). Selection 
effects and differences between the Palomar BQS 
and the Sloan Digital Sky Survey (SDSS) sam- 
ples have been extensively studied by Jester et al. 
(2005). The completeness of the BQS is reported 
to be smaller than that of other surveys such as 
the SDSS and the Hamburg-ESO Quasar Survey 
(Wisotzki et al. 2000). Due to its i?-band mag- 
nitude cutoff, rag < 16, the BQS sample selects 
bright local AGNs, which are likely to be presently 
accreting at high rates (see § 2.1). 

A sample compiled in an opposite manner, i.e. 
with magnitudes fainter than ~16 mags, will prob- 
ably include AGNs that accrete at lower rates than 
PG QSOs. Such a sample is the U-band magni- 
tude my > 15 radio-loud (RL) or radio-quiet (RQ) 
QSO population at 0.138 < z < 0.258 of Dunlop et 
al. (2003). These authors selected objects fainter 
than 15 mags in the V-band, to match the opti- 
cal properties of their entire sample to those of 
their RL subsample (Dunlop et al. 1993). These 
RL QSOs and their optical RQ counterparts have, 
indeed, supcrmassive black holes of 10 8 -10 9 Mq 
that currently accrete at low rates (~0.055 of the 
Eddington value) and are typically located in mas- 
sive elliptical hosts (of mean bulge mass 5.6 x 10 11 
Mq). Tacconietal. (2002), VeiUeux et al. (2002), 
and Dunlop et al. (2003) showed that ULIRGs 
and giant-host QSOs have a very small overlap 



in their host photometric properties and extents, 
which indicates that their origin may be sought in 
different mechanisms. We find that the mean dy- 
namical properties of PG QSOs are intermediate 
to those of these two populations, but significantly 
closer to those of ULIRGs. 

A similar conclusion is derived from the com- 
parison to SDSS radio-loud sources. Best et al. 
(2005) have found a correlation between radio- 
activity and the size of the host galaxy and black 
hole mass; the bulk of radio emission mainly orig- 
inates from the sources with the most massive 
black holes (log[MBH~] ^ 8.2) and the most mas- 
sive bulges in their sample. The latter seem to 
have host-galaxy properties closer to those of the 
Dunlop et al. (2003) sample than those of the PG 
QSOs in this study. 

A possible scenario for the local formation of 5 x 
10 8 -10 9 Mq BH QSOs is by the rcignition of a 10 7 - 
10 8 Mq black hole by the delivery of new gas due 
to tidal interactions. In the case of local galactic 
mergers this delivery could happen, for example, 
via elliptical-spiral or elliptical-elliptical mergers, 
implying that a sequence of merging events would 
be required for the build-up of the black hole. A 
more likely possibility is that these QSOs formed 
at earlier epochs, similarly to some giant Es (see 
Paper II and references therein). 

6. Conclusions 

We have acquired ff-band long-slit spectra of 
12 QSOs, mainly drawn from the PG catalog, to 
study the dynamical properties of their hosts and 
investigate whether they have an origin analogous 
to that of ULIRGs. The compiled sample is rep- 
resentative of local IR-detected PG QSOs and ap- 
propriate for investigating the triggering of QSO 
activity by mergers. We find that: 

1. The long- integration, excellent-seeing-condi- 
tion spectroscopic data obtained from the 
VLT have enabled us to directly extract 
host-galaxy velocity dispersion and dynamic 
properties of local QSOs in the H band. The 
individual stellar velocity dispersions vary 
from 156 to 237 kms -1 . The mean a value 
of the sample is 186 kins" 1 with a standard 
deviation and error of 24 and 7 kms -1 . The 
host galaxies have bulges of ~1.5 m* stellar 
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mass. The dynamically-determined black- 
hole masses of the PG QSOs are of order 
~ 10 s Mq, and they accrete at Eddington 
efficiencies of order ~ 10 _1 . 

2. The measurement of the host velocity dis- 
persion enables us to place 4 reverberation- 
mapped PG QSOs on the AGN M B h - o 
relation, acknowledging the large error bars 
and possible systematics of the measure- 
ment. Most are located above the AGN best- 
fit functions. When using indirect, single- 
epoch, Mbh measurements for 10 QSOs, the 
mean offset from the best fit to the AGN 
data is rather small. The high- and low- 
mass ends of the relation seem consistent 
with each other, implying that this relation 
probably applies to AGNs with a variety of 
Mbh sizes. 

3. The differences in the host-galaxy and black- 
hole properties of our entire sample and the 
IR-excess subsample are insignificant. The 
position of PG QSOs on the fundamental 
plane of early-type galaxies lies between the 
loci of moderate-mass and giant Es. It co- 
incides with the region where the most mas- 
sive ULIRGs are typically located, indicat- 
ing that a subset of PG QSOs and ULIRGs 
have similar dynamical properties. 

4. In some sources, the starburst and strong 
AGN-accretion phase is unambiguously trig- 
gered by a merger event. A scenario where 
some ULIRGs may undergo a QSO phase 
as the merger evolves is therefore plausible. 
The same applies for a scenario which as- 
sumes that some ULIRGs and QSOs may 
have similar, but not identical, progenitors. 
Since there is no evidence that uniquely fa- 
vors one of these scenarios, both may hold 
across the PG QSO population. Still, PG 
QSOs constitute only a fraction of the lo- 
cal QSO population, and sources of Mbh 
5 x 10 8 -10 9 Mq seem to have a different ori- 
gin. 
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Table 1 
QSO List 



Galaxy 


ISA a 


Doc b 




m B 




slit P.A. c 


integration 


J? d 
"off 


< Crff > 




(2000) 


(2000) 




mag 




(°) 


(mins) 


(kpc) 


(mag arc-sec 


PG 0007+106 


00:10:31.0 


+ 10:58:30 


0.0893 


16.1 


f 


139,49 


160,160 


3.66 


17.1 ' 


PG 0050+124 


00:53:34.9 


+ 12:41:36 


0.0611 


14.4 


f 


134,44 


80,80 


1.97 


15.1 ' 


LBQS 0307-0101 


03:10:27.8 


-00:49:51 


0.0804 


16.3 


S 


44,134 


120,120 






PG 1119+120 


11:21:47.1 


+ 11:44:18 


0.0502 


15.4 


f 


89,179 


80,80 


1.39 


15.4 1 


PG 1126-041 


11:29:16.6 


-04:24:08 


0.0600 


15.4 


f 


154,64 


110,105 


4.21 


17.4 1 


PG 1211+143 


12:14:17.7 


+ 14:03:13 


0.0809 


14.4 


h 


0,89 


130,130 


2.57 


16.6 J 


PG 1229+204 


12:32:03.6 


+ 20:09:29 


0.0603 


15.5 


h 


29,119 


140,110 


5.12 


17.6 ' 


PG 1404+226 


14:06:21.8 


+ 22:23:46 


0.0980 


16.5 


h 


45,135 


120,120 


7.15 


18.5 * 


PG 1426+015 


14:29:06.6 


+ 01:17:06 


0.0865 


15.7 


h 


60,-30 


180,120 


6.49 


17.3 > 


PG 1617+175 


16:20:11.3 


+ 17:24:28 


0.1124 


15.5 


h 


-1,89 


120,110 


2.43 


16.5 J 


PG 2130+099 


21:32:27.8 


+ 10:08:19 


0.0630 


14.7 


h 


-41,49 


150,70 


4.44 


17.6 ' 


PG 2214+139 


22:17:12.2 


+ 14:14:21 


0.0658 


15.0 


f 


-1,89 


120,120 


3.97 


17.4 J 



Units of right ascension arc hours, minutes and seconds, and units of declination are degrees, arcminutes, and arcseconds. 
c Slit position angles and respective integration times per slit. 

d The effective radii arc from McLcod & Riekc (1994a) and Veilleux ct al. (2006). 
e The H — band mean surface brightness within the effective radius. 
f Apparcnt magnitudes from Schmidt & Green (1983). 

SValue taken from Chaffee et al. (1991). Sloan Digital Sky Survey data indicate a G— and an R— band magnitude of 15.92 
and 15.77 respectively. An interpolation between these two values yields an expected B-band magnitude somewhat smaller than 
that of Chaffee et al. (1991). 

^The (epoch-averaged) B-band magnitude is from Kaspi et al. (2000). 

1 The K"-band host surface brightness is computed from NICMOS H-band magnitudes {Veilleux et al. 2006) using an H — K 
color index of 0.3 mag (Jarrett ct al. 2003). 

J The it-band host surface brightness is computed from NICMOS if-band magnitudes (McLcod & Rieke 1994a) using an H-K 
color index of 0.3 mag (Jarrett et al. 2003). 
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Table 2 

Source IR fluxes and luminosities 



Galaxy 


log(Ai A [51OOA]/L ) a 


/(12H b 


/(25„m) b 


/(60 M m) b 


/(100 M m) b 


log(i IR /L ) c 


- L IR./ i bbb d 




m.Ty 


mJy 


mJy 


mJy 




PG 0007+106 


11.23 


83° 


163 


192° 


221 f 


11.34 


0.38 


PG 0050+124 


11.21 


549 


1097 


2293 


2959 


11.93 


1.77 


LBQS 0307—0101 
















PG 1119+120 


10.55 


120 


280 


546 


746 


11.13 


0.72 


PG 1126-041 


10.80 


104 


309 


669 


1172 


11.34 


0.69 


PG 1211 + 143 


11.17 


166° 


331° 


412° 


689 


11.58 


0.41 


PG 1229+204 


10.50 


117 


230 c 


202° 


317 f 


11.13 


0.66 


PG 1404+226 


10.80 


26 f 


62 f 


<154 


123 f 


<11.07 h 




PG 1426+015 


11.14 


130° 


221 c 


318 


350 f 


11.50 


0.95 


PG 1617+175 


10.90 


65K 


7lS 


102K 


544 s 


11.43 


0.40 


PG 2130+099 


10.88 


186° 


357° 


480° 


485 f 


11.39 


0.59 


PG 2214 + 139 


11.08 


61 


95 


337 


<282 


<11.04 h 


0.32 



a Optical luminosities taken from Peterson et al. (2004) or Vcstcrgaard & Peterson (2006). 
^All data are from Sanders et al. (1989) unless otherwise noted. 

c Infrarcd luminosities calculated using the Sanders & Mirabel (1996) prescription from the MIR fluxes. 

d Ratio of integrated IR. over big blue bump luminosity (sec text for definition) arc derived from Surace et al. (2001) and Guyon et al. 
(2006). 

e Mean of the Sanders et al. (1989) and Haas et al. (2003) flux measurement. 
f Infrared Space Observatory flux from Haas et al. (2003). 
Slnfrarcd Astronomical Satellite flux from Moshir et al. (1990). 

h ower limit for 1o S (Lj R /Lq) of PG 1404 + 266 is 10.89. when assigning /(60^m) = 0. Lower limit for log(.L /i? /LQ ) of PG 2214 + 139 is 
10.99, when assigning /(100^m) = 0. 
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Table 3 

Host dispersions, black hole masses and Eddington efficiencies 



Source 








apert . a 


template 


«BH b 


M BH (Bdd.) " 


"JEdd d 


M BH (rev.) e 


M 


BH< SE -) 


(PG) 


(kn 






(") 


star 


(«0) 


(M ) 




(Mg) 




(Mq) 


PG 0007+106 


201 


(± 


61) 


0.22-1.10 


M0III 


1.38 X 10 8 


4.05 X 10 7 


0.29 




5 


35 X 10 8 


PG 0050 + 124 


188 


(± 


36) 


0.22-0.96 


KOI 


1.05 X 10 8 


3.85 X 10 7 


0.37 




2 


76 X 10 7 


LBQS 0307-0101 


207 


( = 


49) 


0.22-0.96 


M0III 


1.55 X 10 8 












PG 1119+120 


162 


( = 


28) 


0.22-0.96 


MOIII 


5.79 X 10 7 


8.39 X 10 6 


0.14 




2 


95 X 10 7 


PG 1126-041 


194 


(± 


29) 


0.22-1.10 


KOI+MOIII 


1.19 X 10 8 


1.50 X 10 7 


0.13 




5 


61 X 10 7 


PG 1211 + 143 














3.48 X 10 7 




9.14 X 10 7 






PG 1229 + 204 


162 


(± 


32) 


0.22-1.25 


KOI+MOIII 


5.79 X 10 7 


7.44 X 10 6 


0.13 


7.32 X 10 7 


1 


38 X 10 8 


PG 1404 + 226 


237 


( = 


52) 


0.22-0.96 


KOI 


2.67 X 10 8 


1.48 X 10 7 


0.06 




7 


74 X 10 6 


PG 1426 + 015 s 


185 


( = 


67) 


0.22-1.10 


MOIII 


9.87 X 10 7 


3.25 X 10 7 


0.33 


1.30 X 10 9 


1 


16 X 10 9 


PG 1617+175 


183 


( = 


47) 


0.37-1.40 


MOIII 


9.45 X 10 7 


1.87 X 10 7 


0.20 


5.94 X 10 8 


6 


76 X 10 8 


PG 2130+099 


172 


( = 


46) 


0.22-0.96 


KOI + MOIII 


7.36 X 10 7 


1.79 X 10 7 


0.24 


4.57 X 10 8 


8 


05 X 10 7 


PG 2214+139 


156 


( = 


18) 


0.37- 1.40 


KOI + MOIII 


4.97 X 10 7 


2.84 X 10 7 


0.57 




3 


56 X 10 8 



f 



a The single- sided aperture within which a was measured. 

"Dynamical black hole masses estimated from their relation to the bulge dispersion (Tremaine et al. 2002). 

c Eddington black hole mass, calculated by attributing all of L bo] to the AGN (L^rf^)- 

d Ratio of Eddington over dynamical black hole mass. 

c Reverberation black hole masses from Peterson et al. (2004). 

^Single-epoch, virial black hole masses from Vestergaard & Peterson (2006) with optical luminosities from Boroson Sz Green (1992). 

g The velocity dispersion of the SW (stellar-light-dominated) component of this interacting system is 154 (± 27) kms" 1 : it is derived 
using the MOIII giant. 
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Fig. 1.— Left panel: The optical vs. IR luminosity of PG QSOs. Right panel: The MIR color-index 
diagram used to identify possible transition objects between ULIRGs and QSOs (see Lipari 1994; Canalizo 
& Stockton 2001). The empirically determined position of ULIRGs on this diagram is indicated by a dotted 
box (Canalizo & Stockton 2001). The power-law and black-body lines roughly separate between AGN- and 
starburst- dominated sources. 

In both panels, the QSOs in our sample are plotted as filled circles. For two QSOs in our sample with upper 
limits on their MIR fluxes, arrows indicate how their position would change with a decrease of the respective 
flux. All other Boroson & Green (1992) PG QSOs with exact MIR flux measurements are plotted as open 
circles. This diagram indicates that in terms of both optical and MIR/FIR properties, the QSOs in our 
sample are representative of the local PG population. 



HD 179323 




1.56 1.58 1.60 1.62 
wavelength (/im) 



Fig. 2. — The ii-band spectrum of the KOIab supcrgiant HD 179323 that is used as a template for the 
extraction of the stellar kinematics of some QSO hosts. 
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Fig. 3. — The ii-band spectra of the QSOs in this study. The selected stellar template for each source 
is overplotted as a solid line after being convolved with the Gaussian that represents the LOS broadening 
function of the source. All spectra are shifted to the rest frame. The emission lines occasionally seen at 
1.611 /Ltm and 1.644 jim correspond to Br 13 and [Fell] respectively. 
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Fig. 3. — continued. 
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Fig. 4. — Left panel: The acquisition image of PG 1426+015 clearly shows that it is an interacting pair of 
4.4 kpc nuclear separation. The QSO corresponds to the NE nucleus. 

Right panel: The spectrum of PG 1426+015 for a slit passing through both nuclei. The x and y axes 
have wavelength and spatial dimensions respectively. This image shows the difference between two 10-min 
integrations at different detector positions; this technique is used to remove the sky lines. The number 
counts of the QSO spectrum are one order of magnitude greater than those of the faint nucleus. The velocity 
dispersions of the NE and the SW nucleus equal 185 (± 67) and 154 (± 27) kms -1 respectively. 
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Fig. 5. — The JT-band fundamental plane of early-type galaxies. In all panels, giant boxy and moderate-mass 
disky Es (squares and circles respectively) are from Bender et al. (1992) and Faber et al. (1997). More 
(cluster) Es (open circles) are from Pahre (1999). The QSOs of this study are plotted as stars and those of 
Dunlop et al. (2003) as asterisks. Upper left panel: The <7-i? e ff projection of the plane. Upper right panel: 
The Roff-^os projection of the plane. Middle panel: The 3-d view of the plane. Lower panel: The plane 
projected as in Pahre (1999) for viewing clarity. 
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100 

a (km/s) 



Fig. 6. — The Mbh — cr relation for AGNs with direct Mbh measurements. In this Figure, we plot the virial 
product Mbh/(/) of reverberation experiments (Peterson et al. 2004) versus the host velocity dispersion to 
avoid uncertainties in the mean value of the factor f. The Seyfert 1 AGNs of Onken ct al. (2004) and Nelson 
ct al. (2004) are shown as circles. A few low-luminosity AGNs with upper limits on their Mbh values are 
indicated by arrows pointing down. The QSO datapoints are shown as stars. The solid and the dashed lines 
correspond to straight lines with slopes identical to those of the quiescent-galaxy relations (Tremaine et al. 
2002 and Ferrarese 2002 respectively) that fit the AGN datapoints (Onken et al. 2004). 




Fig. 7. — The AGN Mbh — o relation studied at its high- and low- mass end using indirect, single-epoch, 
Mbh estimates for local PG QSOs (Vestergaard & Peterson 2006) and low-mass AGNs (of Mbh ~ 10 6 ; 
Barth et al. 2005). As in Fig. 6, the value (/) = 5.5 is divided out of all black-hole estimates. The symbols 
used for the Seyfert 1 AGNs of Onken et al. (2004) and Nelson et al. (2004) and the PG QSOs of this study 
are identical to those used in Fig. 6. The low-mass AGNs are plotted as triangles. The solid and the dashed 
lines correspond to the Tremaine et al. (2002) and Ferrarese (2002) fits scaled to the AGN datapoints. 
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Fig. 8. — The X-band fundamental plane of merger remnants. The panels are identical to those of Fig. 5. 
In all panels, ULIRG, LIRG, and other (visually-selected) remnants are plotted as triangles, circles, and 
squares respectively. The merger remnant data are from Genzel et al. (2001), Tacconi et al. (2002), Paper 
II, Rothberg & Joseph (2006), Hinz & Rieke (2006), Shier & Fischer (1998), and James et al. (1999). The 
QSOs of this study are plotted as stars and those of Dunlop et al. (2003) as asterisks. 
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